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Abstract 

We consider thermal n-point Green functions in the framework of quan- 
tum field theory at finite temperature. We show how analytic continuations 
from imaginary to real energies relate these functions originally defined in the 
imaginary-time formalism to retarded and advanced real-time ones. The de- 
scribed method is valid to all orders of perturbation theory. It has the further 
advantage that it is independent of approximations often applied in actual 

finite-order calculations. 
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I. INTRODUCTION 



In quantum field theory at finite temperature and density two convenient formalisms 
that enable the use of conventional Feynman rules in momentum space are applied; the 
imaginary-time formalism (ITF) and the real-time formalism (RTF) [0,0. (For books 
and reviews on thermal field theory, see, e.g. The former one is particularly suited 

for the evaluation of the static or thermodynamic quantities of finite-temperature systems, 
while the latter is preferred for the evaluation of time- dependent quantities. 

Studies of the two-point functions have a long history ||T]-0. The relationship between 
their representations — the one obtained from ITF and the other from RTF — is well known 
1^. On the other hand the interest in the three- and n{> 4)-point functions started rather 
recently [^-[14 . 

An apparent difference |jlO| between the results for three-point functions obtained from 
ITF and from RTF has posed the question: through analytic continuations of the Green 
functions evaluated in ITF, what kind of (combinations of) thermal Green functions in RTF 
are obtained? 

Several papers have been devoted to this issue. It has been shown, either to one-loop 
order or in all orders of perturbation theory, that the three-point function in ITF, when 
analytically continued to real energies, becomes the retarded or advanced three-point Green 
function ||TIl , |l2| . Recently, the n (> 4)-point functions have been investigated jTSlJl^ and it 
is shown that different analytic continuations of the ITF result yield different RTF Green 
functions, including the retarded and advanced Green functions. 

The purpose of this paper is to show that the most straightforward analytic continua- 
tion of the non-amputated ra-point Green function in ITF leads to the retarded or advanced 
n-point Green function. 

In Sec. II, we introduce the n-point thermal Green functions with time arguments on 
a — to a large extent arbitrary — contour in the complex time plane. We then formulate 
the problem of analytic continuation of the Green function in ITF from imaginary- to real- 
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energies. For the purpose of illustrating our procedure, we discuss in Sec. Ill the analytic 
continuation of the two-point function. In Sec. IV, we carry out the analytic continuation of 
the n-point Green functions evaluated in ITF, and obtain the retarded- and advanced- Green 
functions. Sec. V is devoted to conclusions. 



II. PRELIMINARIES 



Throughout this paper we consider a real scalar field 4>{x). Generalizations to other 
kind of fields are straightforward. The thermal Green functions are defined as the statistical 
average of a product of Heisenberg fields, 

G{{t})^G{h,t2,---,tn) 



= Tr 



= {T,{<t>{t,)<t>{t2)---(t>{tn))) , 



(1) 



where (3 — T"^ is the inverse temperature and H is the Hamiltonian of the system such that 
= e*^*0(O)e~*^*. In Eq. (1) and in the following we suppress explicit reference to the 
space variables. The arguments ii, t2, - ■ ■ ,tn lie on the contour C running from to to to — ij3 
in the complex time plane. The symbol Tc in (1) is the time-ordering operator along this 
path C. That is, it prescribes that the operators it is applied to be arranged in the order 
in which their arguments lie along C, with those nearest to the beginning at to to the right, 
and those nearest to the end to — iP to the left. 

To perform the thermal trace (1), we insert a complete set of states by choosing the 
Hamiltonian eigenstate basis: 



G{{t})Tre-^" = ^ 



Pn 
X 



n— 1 



n^c(tj-t^) 



^ exp{-iEi^ {tn -tj-ip)) 



j=2 



{ln+l = ll), (2) 



with 6c the contour step function. The summation here is carried out over all permutations 
of n numbers p„ 
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Following common practice, we assume that the convergence of the above trace sum (2) 
is controlled by the exponential factors, so that it converges if, for every pair of ti and tj 
such that edti - tj) = 1, ^{ti - tj) < 0, and ^{ts - ti) < (3 with t^ (ti) the "smallest" 
("largest") time. This condition guarantees the existence of G{{t}) (Eq. (2)) as an analytic 
function of {t} = {ti, t2, - ■ ■ , tn}- The limit of an analytic function on the boundaries of its 
domain of definition, where it is still continuous, is a generalized function. This implies that 
the thermal Green function G{{t}) is well defined for S>(tj — tj) < when 9c{ti — tj) — 1, 
and '^{ts — ti) < (3. This imposes the restriction on C that as a point moves along C from to 
to ^0 ~ its imaginary part is nonincreasing. Then, an analytic continuation of G{{t}) can 
be done by deforming the contour C with the end points to and to — iP held fixed, keeping 
in mind the "nonincreasing" condition for the imaginary part of the points on C. 

Among the above class of paths, we consider first a special path Cj depicted in Fig. 1, 
which defines the ITF. We note that we can choose any value for 3? to because of the property 
of time-translation invariance of (1). Next we evaluate the Fourier component of (1), 



It is to be noted that real (discrete) utj's here are what we call imaginary energies. By using 
once more the time translation invariance we rewrite (4) as 




,=1 \JCr / 
ujj = 27rlj//3, = 0,±l,±2,---,±oo . 



(4) 



n 



G{{u;}) = -i(3 SiY^ Uj ; 0) G{u2, u^,---, 



(5a) 



G{llI2, UJ3, ■ ■ ■ 




(5b) 
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Here the integrations are performed along the imaginary-time axis. In (5a) S{- ■ ■ ; ■ ■ ■) 
denotes the Kronecker 6 symbol. In order to obtain (5) the KMS condition [|I5|, which 
represents the invariance of the trace under the following cyclic permutation (cf. (1)), 

■ ■ ■ (f){tn-l)(j)itn)) = (0(t„ - ■ ■ ■ (6) 

is used. In ITF one calculates G or G as defined in (5). In the following we focus our 
attention on how to continue (5) from imaginary to real energies. 

III. TWO-POINT THERMAL GREEN FUNCTION 

Although the relation between ITF and RTF is well understood ^ for two-point thermal 
Green functions, we start with the thermal two-point Green function in ITF for the purpose 
of illustrating our procedure: 

= [ dtj dt2e-^^^'^+'^'''\TcMtimh))) . (7) 

JCi JCi 

As explained in Sec. II, the integrand of (7) is an analytic function of ti and t2 in the 
strip, '^to — P < '^tj < 53 to (j = 1,2) with ti ^ t2, and we may deform the contour Cj 
keeping the property as mentioned above after (3). In this way we may obtain the contour 
Cr = Ci © 6*2 © C3 © C4 as depicted in Fig. 2: the path goes down from to to t/ = 3? to, 
continues from tj to tp along the real axis, returns back to tj along the real axis, and ends 
at to — ip. 

Because of time translation invariance G{uji,uj2) is nonvanishing only for cui + ct;2 = 
(cf.(5a)), and the integrand of (7) — on the path Cr — is a function of t2 — ti. This, together 
with the above mentioned analyticity property of G(ti,t2), allows us to evaluate 0(001,002) 
by fixing ti on any point on the contour Cr = Ci © C2 © C3 © C4. We first fix ti on Ci, to 
obtain 



C{ooi,o02) = -il36{ooi; -002) 



rft2e--^(*^-*^)(0(ti)0(t2)) 



+ / rft2e-"^(*^-*^)(Tc,(0(ti)0(t2)))+ / cit2e-"^(*^-*^)(0(t2)0(ti)) (8a) 

= -zP Siui; -u^) r dh e-'^^(*^-*i) [(T(</)(ti)</)(t2))) - m^mh))] 
Jti 

+ (contributions from C3 and C4) . (8b) 

The symbol T is the ordinary time-ordering symbol. Eq. (8) is now well suited for the 
purpose of an analytic continuation of G{u!i,u!2) to real energies. In order to arrive at RTF 
with real energy p2o{= —pio) we take the limit ti{= 5Rto) —00 and tp +00, taking into 
account that ti{= Kto) niay be chosen arbitrarily. 

We realize (cf. 0) that the term in the square bracket in (8b) may be written as 

(r(0(ti)0(t2))) - (0(t2)0(tl)) 

= ^(tl-t2)([0(tl),0(t2)]) . (9) 

Then, from (8) and (9), we learn that the above limit, ti —00 and tp +00, can be 
taken if we continue to the real energy as follows, 

UJ2 -iiP20 - «e) , (10a) 

or equivalently, 

ui -i(pio + ie) , (10b) 

with e an infinitesimal positive number. 

The time path for G consists of two real-time segments, namely Ci from —00 +00 
and C2 from +00 to —00, and of C3 © C4 (C3 from —00 + i'^ to to —00 and C4 from —00 to 
—oo + iC^to — P)). Then, in evaluating (8) with (10) by perturbative methods, there are four 
kind of bare thermal propagators, Gf^{ti, tj) = 1 — 4), present with ti e Cj. Following 
||16|| , one can show that as far as the thermal Green functions (like (8) with (10)) with their 
time arguments lying on Ci and/or C2 are concerned, the net effects of the contributions 



from C3 and C4 are to modify G^f = 1, 2) to the familiar forms 0] in two-component 
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real-time thermal field theory constructed on the basis of the closed-time path Ci © C2 
ill- 

The analytic continuation of the energy conservation 5-function in (8) is given by the 
prescription 0: -CJ2) 27rz/5~^5(pio + P2o)- 

Thus we obtain the analytically continued G with real energies pio and P20, 

lim G{pio + it, P20 - it) = G++ - G+- (11a) 
= 2nS{pio + P20) 

/OD 
dte'^P"^'e{-t){[(f){0),(f){t)]) , (lib) 
-00 

where we follow e.g. p and introduce 

/OO nOO 
dtJ rft2e^^^^°*^+P-*^)(A„^(ti,t2)), (12) 
-00 J —00 

with the definitions 

A++(ti, ta) = 0{ti - t2)(j){ti)(j){t2) + e(t2 - h)Ht2)Hti) 

A+_{t^,t2) = (f){t2)(j){ti) 

A^4h,t2) = e{h - ^2)0(^2)0(^1) + 0{t2 - ti)0(ti)0(t2) 

A_+{tut2) = (f){h)(j){t2) . (13) 

The functions Gais (a,/5 = +, — ) in (12) denote Green functions in the two-component real- 
time thermal field theory, applying the closed time-path formalism introduced by Schwinger 
and Keldysh [^. It is to be noted, in passing, that the A^/j's are not independent since they 
satisfy the relation 0, 

A++{h,t2) + A__{tut2) 

-v4+„(ti,t2)-v4_+(ti,t2) = 0. (14) 
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It is important to note that the continued function G in (lib) is identical with the retarded 
Green function. In this way the "physical" representation as discussed in is established, 
where this function is denoted by G = iG2i{pw,P2o), however, the primary quantities that 
are calculated directly in real-time thermal field theory are Gap in (12). 

When we fix instead of ti, on Ci, and repeating similar steps as above we obtain 
the retarded Green function with respect to 0(^2), i-e. the advanced Green function with 
respect to 0(ti). Fixing the time variable as ti G G2 leads to the same result as above, (11), 
and likewise for the choice t2 G C2. Fixing ti on C3 or C4, ti G C3 © C4, is not suitable 
for analytic continuations under consideration and we recover the original formula (7); this 
holds likewise for the choice ^2 ^ C3 © C4. 

Finally we may deform the contour C/ in (7) to the one that is mirror symmetric to the 
one of Fig. 2. Changing the time variable tj as tj = ^tj + i'^tj — 3?tj + i'^tj, we get 
back the time path Ci © C2 © C3 © C4 of Fig. 2. Then, fixing ti on the upper path on the 
real axis in the complex time-plane, we deduce 

lim G(pio - ie, P20 + ie) = -G++ + G_+ (15a) 

e— >+0 

= -2nS{pio+p2o) 

/oo 
c?te<^-+-)*^(t)([0(O),0(i)]), (15b) 
-00 

i.e. the advanced Green function. 

IV. THERMAL A^-POINT GREEN FUNCTION 

We proceed in a similar manner as in Sec. III. Starting with the thermal n-point function 
in ITF, Eq. (4), we deform the contour Gj to Gr as depicted in Fig. 2. Under the constraint 
= (cf. (5a)), we evaluate (4) with Gr for Gj by fixing ti on Gi, ti G Ci, which is 
suitable for continuation to real energies. In place of Eq. (8) we now have 

Giico}) = -tpSiY.^,;0) ft f / dt,e-^^'^-'A 

j j—2 \"'CiffiC2 / 

8 



X (Tc,eC2(0(tl)---0(tn-l)0(tn))) + 



(ti e Ci) , 



(16a) 



1 1=5 ^"'^i^ 



X 



E (-)'(rC,eC. (0+(tl)0a,(t2) ■ ■ •0a„(tn))) 



+ 



(16b) 



where 



1 



E(i- 



a 



3J ' 



(17) 



and the notation of is used: for n > 2 it becomes more transparent to make the subscripts 
+, — exphcit, which indicate that the times ti assume values on either the "positive" time 
path Ci or on the "negative" one C2; e.g. this imphes for all + subscripts time ordering, 
whereas for all — subscripts anti-time ordering. This generalizes the expressions in (13) to 
the n > 3 cases (for n=3 explicit expressions are given in 0). In Eq. (16) the dots indicate 
the contributions when some of the t's among t2,---,tn are on C^^Ci, this portion is treated 
according to the arguments given in Sec. Ill ||16 |. 

Next we transform to the "physical" representation. By applying the procedure described 
in 0, we express the term in the square bracket in (16b) as, 



E^(l,2,---,n)([---[0(ti),0(%)],0(%)], 



0(k)]) 



(18) 



Pn-l 



where 9 is the multi-step function defined by ^^(1, 2, ■ ■ ■ , n) = ^^(1, 2) 9(2, 3) • ■ ■ 9{n — 1, n) 
with ^(1,2) = 9(ti — 12). Eq. (18) is the generalization of the two-point function case Eq. 
(9). The summation here is carried out over all permutations oi n — 1 numbers Pn-i 

\ 



2 3 ■■■ n 
2 3 ■ ■• n 



(19) 



Eq. (18) tells us that ti is the largest time, therefore we are allowed to take the limit 
ti{= 3?to) ~oo and tp — » +00 if we continue Uj in (16b) to real energies as 



Uj -iiPjo - ie) , j = 2, 3, 



(20a) 



and 



oji -i{pio + i{n - l)e) . (20b) 
Then we arrive at the Green function in RTF, 

^lim^ G{pio + i{n - l)e, {pjo - ie; j = 2, 3, ■ ■ ■ , n}) = (-)'G'+,02,-,«n (21a) 

X ^ e{l,% ■■■,n){[--- [0(ti), 0(%)], </)(%)], • ■ ■ ], 0(t^)]) (21b) 

Pn-l 

= r^^Gaii.-.i (pio, ■ ■ ■ , Pno) ■ (21c) 

where s is given in (17). Thus we have derived the Green function (21c) in the "physical" 
representation, denoted by G2ii - i lH, which is the thermal n-point retarded function as seen 
in (21b). In (21a), Gj^^a2,---,a„ is a Green function in RTF and it is defined @] analogously 
to (12) and (13), 



/ roc \ 



X (Tc,eC2 (0o,(tl)0a,(t2) ■ ■ ■0a„(tn))) • (22) 

As in the two-point function case, (14), there is one identity; 

X! Gai,a2,-,an = (23) 

where 

^' = ^E(l-«.)- (24) 

^ i=i 

It is worth mentioning that the primary quantities that are evaluated in real-time thermal 
field theory are ,02, defined in (22), through which the retarded Green function G in 
(21b) is obtained. 
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We have developed the derivation by fixing ti on Ci. Of course, we may proceed in a 
similar manner by fixing other tj (2 < j < n) on Ci: n — 1 different retarded Green functions 
are the result. 

In case Eq.(4) is evaluated with Cr for Cj by fixing ti on C2, we obtain the same result 
as above, (21); for ti G C3 © C4 we recover the original formula (4). 

When we deform the contour Cj in (4) to the one that is mirror symmetric to the one of 
Fig. 2, and fixing ti on the upper path on the real axis in the complex-time plane, i.e. the 
counterpart of Ci in Fig. 2, we derive the advanced Green function, 

^lim^G'(pio - i{n - l)e, {pjo + ie; j = 2, 3, ■ ■ ■ , n}) 

n " / /-oo \ 

= -2Ti5{Y,Pjo) lim n / rftje^(f^o+*^)(*^-*iM 

X ^(^, ■ ■ ■ , 2, 1)([ - ■ ■ 0(%)], 0(%)], ■ ■ ■], . (25) 

Pn-l 

V. CONCLUSIONS 

In this paper we have addressed the question: what kind of thermal functions in RTF 
emerge by analytic continuations of the ra-point thermal Green functions defined in ITF. 
This amounts to perform analytic continuations in the energies of the external legs from the 
discrete imaginary values to real continuous ones. 

The thermal Green functions are defined on a path in the complex time plane, a path 
which is to a large extent arbitrary. On the basis of this observation, we have carried out 
the above mentioned continuations in the most straightforward and familiar manner by 
deforming the contour, starting from the one that defines ITF to the one defining RTF. In 
this way, we show that ITF n-point Green functions become retarded or advanced thermal 
Green functions. 

The results obtained in this paper, being exact and valid independent of the approxi- 



mation used in actual mainly perturbative calculations, are partly discussed in |13], and we 
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expect that our approach gives additional information on the underlying relations between 
ITF and RTF. 



12 



REFERENCES 
T. Matsubara, Prog. Theor. Phys. 14, 351 (1955). 

Y. Takahashi and H. Umezawa, Collet. Phenom. 2, 55 (1975); A. J. Niemi and G. W. 
Semenoff, Ann. Phys. (N.Y.) 152, 105 (1984); Nucl Phys. B230, 181 (1984). 

J. Schwinger, J. Math. Phys. 2, 407 (1961); L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 
1515 (1964) [Sov. Phys. JETP 20, 1018 (1965)]. See also, R. A. Craig, J. Math. Phys. 
9, 605 (1968). 

A. A. Abrikosov, L. P. Gor'kov, and I. E. Dzyaloshinskii, Methods of Quantum Field 
Theory in Statistical Physics (Pergamon, Oxford, 1965); H. Umezawa, H. Matsumoto, 
and M. Tachiki, Thermo Field Dynamics and Condensed States (North- Holland, Ams- 
terdam, 1982); J. I. Kapusta, Finite-temperature Field Theory (Cambridge University 
Press, Cambridge, England, 1989); 

R. Mills, Propagators for Many-Particle Systems (Gordon and Breach, New York, 1969). 
K.-C. Chou, Z.-B. Su, B.-L. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 
N. P. Landsman and Ch. G. van Weert, Phys. Rep. 145, 141 (1987). 
G. Baym and N. D. Mcrmin, J. Math. Phys. 2, 232 (1961). 

E. Braaten and R. D. Pisarski, Nucl. Phys. B337, 569 (1990); J. Prenkel and J. C. 
Taylor, ibzd. 334, 199 (1990). 

R. Baier, B. Pire, and D. Schiff, Phys. Lett. B238, 367 (1990); Z. Phys. C51, 581 
(1991); H. Nakkagawa, A. Niegawa, and H. Yokota, Phys. Lett. B244, 63 (1990); 273, 
549(E) (1991); N. Ashida, A. Niegawa, H. Nakkagawa, and H. Yokota, Phys. Rev. D44, 
473 (1991); 45, 1432(E) (1992). Earlier work is quoted in these papers. 

T. S. Evans, Phys. Lett. B249, 286 (1990); 252, 108 (1990). 

R. Kobes, Phys. Rev. D42, 562 (1990); 43, 1269 (1991); Phys. Rev. Lett. 67, 1384 

13 



(1991); P. Aurenche, E. Petitgirard, and T. del Rio Gaztelurrutia, Phys. Lett. B297, 
337 (1992). 

T. S. Evans, Nucl. Phys. B374, 340 (1992). 

P. Aurenche and T. Becherrawy, Nucl. Phys. B379, 259 (1992); J. C. Taylor, Phys. 
Rev. D47, 725 (1993). See also, M. A. van Eijck and Ch. G. van Weert, Phys. Lett. 
B277, 305 (1992); F. Guerin, Institut Non Lineaire de Nice Report No. INLN 1993/14. 

R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957); P. C. Martin and J. Schwinger, Phys. Rev. 
115, 1342 (1959). 

A. Niegawa, Phys. Rev. D40, 1199 (1989). 

The modifications are as follows 0,0: gS? ^ G++, gS? ^ G'+_, Cg^ ^ G^^ and 
G21 Here G's are given in (12) and (13) below. 



14 



FIGURE CAPTIONS 
FIG.l. The contour Cj in the complex time plane, which defines ITF. 
FIG. 2. The contour Cr = C3 © Ci © C2 © C4 in the complex time plane; the segments Ci 

and C2 lie on the real axis, and ti — 3?io- 
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